Universal Control of Decoupled Quantum Systems 
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It is shown that if one can perform a restricted set of fast manipulations on a quantum system, 
one can implement a large class of dynamical evolutions by effectively removing or introducing 
selected Hamiltonians. The procedure can be used to achieve universal noise-tolerant control based 
on purely unitary open-loop transformations of the dynamics. As a result, it is in principle possible 
to perform noise-protected universal quantum computation using no extra space resources. 
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The desire to shape quantum evolution according to 
precisely controlled dynamics is shared by various areas 
of contemporary physics and engineering [pj. A prob- 
lem commonly encountered in the task of controlling the 
dynamical behavior of a quantum system is the need 
for removing unwanted interactions present in the full 
Hamiltonian: historically, one of the first solutions was 
provided by Nuclear Magnetic Resonance spectroscopy, 
where a variety of decoupling techniques have been de- 
veloped to simplify spectra by effectively elinrinating se- 
lected contributions to the nuclear Hamiltonian M . Once 
suppression of a given term is obtained, the correspond- 
ing Hanriltonian may no longer be directly available for 
control. From the perspective of attaining universal dy- 
namical control, this raises the question of devising ways 
for introducing or re-introducing Hamiltonian control 
compatible with the prescribed decoupling action. 

A striking example is the case of an open quantum sys- 
tem, where the coupling to the environment is responsible 
for inducing quantum decoherence and dissipation pro- 
cesses, thereby corrupting the original unitary dynamics 
H . The demand for universal control strategies of noise- 
decoupled quantum systems, able to effectively reject en- 
vironmental noise while still ensuring full control capabil- 
ities, has been heightened tremendously due to the chal- 
lenge of implementing quantum computation B . In spite 
of the beautiful and powerful advancements made in the 
theory of quantum error correction M , fault-tolerant er- 
ror correction || and concatenated coding Q, practical 
exploitation of these results is still seriously constrained 
by the amount of extra space resources required || . 

In this Letter, we address the general problem of open- 
loop controllability of a decoupled quantum system: the 
controller is assumed to apply time-dependent potentials 
without ever measuring the actual state of the system. 
We introduce programming procedures for combining the 
desired control action with the decoupling operations and 
identify the conditions under which universal control over 
the effective decoupled dynamics is retained. As a conse- 
quence, we demonstrate the possibility of achieving noise- 
tolerant control of open quantum systems solely based on 
unitary manipulations which do not require ancillary re- 
sources. 



Decoupling. — We begin by recalling the essential in- 
gredients of decoupling in the language of J9| . Let H be 
the Hamiltonian of a quantum system living in a Hilbert 
space TL. Suppose we have the capability of performing 
instantaneously a certain set of unitaries, meaning that 
the corresponding set of Hamiltonians can be turned on 
for negligible amounts of time r with (ideally) arbitrarily 
large strength. We shall term such impulsive full-power 
control operations as bang-bang (b.b.) controls PJlOf. If 
U is a rotation that can be implemented b.b., it is con- 
ceivable that U^ 1 — W can be realized b.b. as well. 
In this case, the set of realizable b.b. operations is a 
subgroup Gb.b. of the full group U(H) of unitary transfor- 
mations over Ti. A decoupler on Ti operates by iterating 
the system through a cyclic time evolution which judi- 
ciously combines sequences of b.b. operations with free 
evolutions under the natural Hamiltonian H. Accord- 
ingly, a decoupler shall be characterized by a finite group 
Q C Gb.b of b.b. operations (decoupling group) together 
with a known time scale T c (cycle time) determining the 
duration of a single cycle. Full knowledge of the decou- 
pler operations, including the exact sequence of group 
operations gj, j = 0, . . . , \G\ — 1, \G\ = ord(Q), and their 
temporal separation At, may not be available from the 
beginning (black-box decoupler) . 

What does a decoupled evolution look like ? A conve- 
nient picture is provided by average Hamiltonian theory 
plpj. Consider a slicing of a given evolution time T as 
T = NT C = N\G\At, and let U Q (At) = exp(-iHAt) be 
the free propagator. The evolution in the presence of 
the decoupler is uniquely determined by the sequence of 
group transformations {gj} over a single cycle time: 
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U ( T c) = II 9)Uv(At)g 3 = e- iH ^^ , (1) 

j=o 

H e ff denoting the resulting effective Hamiltonian. In 
the ideal limit of arbitrarily fast cycle time T c — ► 0, with 
N — > oo in such a way that NT C = T , H e ff approaches 



H » H eff = 



= W\T,9)Hg 3 ^ g (H) 



(2) 
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Eq. (g) defines a quantum operation on the space 
End(H) of bounded operators on H. Physically, H e ff 
represents the leading contribution to the average Hamil- 
tonian describing the motion of the system under the 
influence of the decoupling field, the j-th term in the 
sum (g) being identical to the so-called toggling frame 
Hamiltonian during the j-th cycle subinterval 0. From 
a geometric point of view, the map Ug introduced in (g) 
can be identified with the projector on the so-called cen- 
tralizer Z(G) in End(H) |Q: 

Z(G) = { O | [O, 9j ] = 0V ft e6} = n e (End(H)) . (3) 

Z(Q) is a subalgebra of End(7Y). We shall denote by 
Zh{Q) the subspace of Hermitian operators belonging to 
Z(Q). Eq. (||) allows a direct interpretation of the de- 
coupler action in terms of symmetry properties: since 
H e ff — Hg(H), the decoupled evolution is symmetrized 
according to the group Q O] . All the components of the 
dynamics generated by H, which are not invariant under 
the group Q, are effectively averaged out. 

The group-theoretical prescription (Q) translates di- 
rectly to pulse control, a control cycle involving a se- 
quence of decoupling pulses Dj = gjg i _ 1 , j = 1, . . . , \Q\, 
separated by delays of free evolution At and fulfilling 
D\g\ . . . D2D1 = 1 by cyclicity. The limit T c — > cannot 
be met exactly in practice. While symmetry properties 
suffice to specify the output of a given decoupler, time 
scales are crucial in determining how good the decoupler 
performs in a realistic scenario where both the pulse du- 
ration r and the cycle time T c are finite. Qualitatively 
speaking, the projection on the centralizer (0) will tend to 
symmetrize interactions whose typical correlation times 
are long on the time scale determined by T c . If r c rep- 
resents the shortest correlation time associated with the 
unwanted interactions, decoupling will be effective under 
the hierarchy of time scales r <C At < T c <C t c Q . Ef- 
fects due to finite pulse duration are expected to be of 
order 0(t/t c ), while cycle time corrections scale with a 
ratio 0(T c /t c ) [91. It is worth stressing that T c deter- 
mines the minimum time scale over which decoupling is 
guaranteed. Thus, the original continuous-time dynamics 
under H is replaced by a stroboscopic time development 
under H e ff with a natural time unit equal to T c . 

For open quantum systems, the above description com- 
prises two coupled subsystems S and B with associated 
Hilbert space H — Tts ® Hb, the S-component repre- 
senting the system of physical interest (e.g., the com- 
putational degrees of freedom of a quantum computer). 
If H = Hs + Hb + Hsb denotes the total Hamiltonian, 
environmental noise is introduced through a set of (trace- 
less) error operators E a in the interaction Hamiltonian, 
Hsb — J2 a Ea <8> B a , B a denoting bath operators. We 
require the linear space £ — span{£" Q } to be finite di- 
mensional. An error space £ will be correctable by any 
decoupler ensuring that Ug(£) = 0. The appropriate 
correlation time r c is related to the memory time of the 
reservoir B pj. Note that this necessarily implies a non- 
Mar kovian error scenario. 



The capability of successfully decoupling the system S 
from the effects of the environment is only a first step to 
achieving noise-tolerant control of S'-degrees of freedom. 
Can we do better? 

Programming.— The presence of a decoupler imposes 
restrictions on the manipulations available to further con- 
trol the system dynamics. Those restrictions are in the 
form of both symmetry and timing constraints. Suppose, 
for instance, that a Hamiltonian A is applied over a time 
interval long compared to T c . Then the evolution gener- 
ated by A may be quenched completely if A $. Zh{Q). If, 
instead, the action corresponding to A is implemented in 
the form of b.b. rotations, it is easy to check that, unless 
the insertion points are chosen carefully, any extra pulse 
interferes with the decoupler operations, thereby spoiling 
noise averaging. This can be avoided, for instance, by 
making sure that the b.b. operations are inserted in be- 
tween cycles only. In general, indiscriminate application 
of control operations will not lead to the expected result. 
The goal of programming is to characterize the degree 
of control attainable compatibly with the desired decou- 
pling action. This is influenced by three factors: avail- 
able knowledge of the error space £; available knowledge 
of the decoupler operations; available control resources. 
In particular, since switching on/off strong interactions 
for short times is difficult in practice, a relevant crite- 
rion will be trying to keep the number of required b.b. 
operations to a minimum. 

Let us first consider the case where the error space £ 
is known and a decoupling group Q exists with a non- 
trivial centralizer, Z(Q) =/= {!}, a situation correspond- 
ing to selective averaging according to ||. By virtue of 
(0) , it is always possible to apply slowly any Hamiltonian 
A 6 Zjj{G) in parallel with the decoupler. However, since 
the unitary evolutions generated by such Hamiltonians 
also lie by construction in the centralizer of Q, universal- 
ity on the full Hilbert space cannot be achieved by purely 
exploiting this kind of weak (strength)/ slow (switching) 
control. Supplementary coding methods are demanded, 
which shall be discussed elsewhere 12 . Instead, we fo- 



cus here on the idea of combining slow control from the 
set Zh(G) with suitable fast manipulations available in 
addition to the decoupling ones. 

A relatively straightforward situation occurs when the 
group Gb.b. of realizable b.b. operations is large enough to 
accommodate two distinct decouplers with known (pos- 
sibly different) cycle times, i.e., Tlg(£) — 0, ILj(£) = 0, 
with Q, Q C Gb.b. ■ Notice that if G and P are respectively 
a decoupling group and a unitary transformation, P ^ G, 
P $ Z(G), then G = P^GP implements a twisted decou- 
pler of the same order provided P^£P — £. Suppose now 
we can apply Hamiltonian A G Zh{G) for a time interval 
ATi = N\T C \ (in parallel with decoupler 1), Hamiltonian 
B <G Zh{G) for a time interval AT2 = ^2T c i (in parallel 
with decoupler 2), etc., (e.g., one could have A — Ug(H), 
B = Wq(H)). Then, by using standard universality re- 
sults |L3J, any U = e L could be created, where L belongs 



to the Lie algebra generated by iA, %B under commuta- 
tion. Accordingly, universal control over the decoupled 
dynamics is achieved whenever this algebra amounts to 
the whole Lie algebra of anti-Hcrmitian operators. 

Even if only a single decoupler Q is available, a similar 
strategy can be mimicked through a simple trick. Sup- 
pose that, in addition to decoupling pulses in G, we can 
perform on the system a b.b. rotation P G Gb.b.- What 
can we do with this capability ? Assuming that we are 
able to synchronize operations with the cycle time, we 
can make the system effectively evolve according to a 
transformed average Hamiltonian. Let AT = NT C be a 
given time window, with decoupled evolution ruled by 
the Hamiltonian (0). Imagine now inserting a pulse P 
immediately before the beginning of AT, followed by a 
pulse P' synchronized with the end of AT. Then evolu- 
tion over AT can be described in terms of a new average 
Hamiltonian H e ff = P j 'Hg(H)P i.e., 



Heff = 



c ^9}Hg^n^H), 



~9,£S 



(4) 



where H = P^HP and a twisted decoupling group 
Q = P^QP has been defined, with associated central- 
izer Z(Q) = P^Z{G)P. Thus, the net effect of the two 
programming pulses P, P^ amounts to implement decou- 
pling according to G, noise averaging being retained since 
both the original decoupling group Q and the error space 
£ are simultaneously rotated. Clearly, one has to ensure 
that P <£ Z(Q) in order to steer the effective Hamiltonian 
out of the original centralizer Z(Q). 

Let A now be, as above, a realizable Hamiltonian in 
the centralizer of Q and let B = P^AP denote its ro- 
tated counterpart. Then, by alternating evolution peri- 
ods according to A and B, the latter being obtained by 
inserting pairs of pulses P, P' with appropriate timing, 
it is possible in principle to obtain any Hamiltonian in 
the algebra generated by A, B under commutation. The 
reasoning is easily extended to the case where a given 
choice of interactions is realizable in Zh(G)- In practice, 
the advantageous feature of this scheme is that by per- 
forming a single extra b.b. operation, a new repertoire of 
Hamiltonians becomes effectively available for slow con- 
trol in the centralizer Z(Q). In the generic case, under 
the conditions given in |13j , any desired unitary transfor- 
mation will be reachable in principle, thereby implying 
complete control of the decoupled evolution. 

Note that knowledge of the exact decoupling sequence 
has not been exploited so far, implying validity of the 
previous schemes even for black-box decouplers. If de- 
tailed information on the decoupler operations is avail- 
able, this knowledge can be used to devise alternate con- 
trol schemes implying less stringent resources. Suppose 
that we want to reintroduce control by some Hamiltonian 
B £ Zh(G). If e lB e Gb.b., one could always, in princi- 
ple, exploit the freedom of inserting such a b.b. pulse at 
the beginning and/or the end of decoupling cycles with- 
out affecting decoupling itself. Actually, it is possible to 



replace the the b.b. requirement with a weaker assump- 
tion, by imagining that the strength of B cannot be made 
arbitrarily large but B can still be turned on and off arbi- 
trarily fast. In other words, let us assume a form of weak 
(strength)/ fast (switching) control whereby Hamiltoni- 
ans can be modulated at the same rate as the b.b. control 
within a cycle. Then to reintroduce control according to 
B over a time interval AT = NT C , it suffices to turn it on 
during the 1-frame subinterval of each decoupling cycle. 
The evolution is ruled by the effective Hamiltonian 



H, 



-// = e(£*}* 



a j eg 



,,,+B) = U g (H) + ^-B, (5) 



which acquires a component along B g" Zh{G)- Strength 
reduction for such a Hamiltonian can be avoided if an 
enlarged set of interactions is amenable of fast switch- 
ing: one just turns on a Hamiltonian Bj = gjBg^ during 
the j-th subinterval of each cycle, the overall effect being 
elimination of the |C/| _1 -factor in front of B. 

Using the above methods, controlled evolutions can be 
designed by both letting the system evolve under the ac- 
tion of the decoupler alone and by incorporating modified 
decoupling cycles to displace the effective Hamiltonian 
out of Z(G) as in (0). The issue of complete control can 
now be addressed by looking at the combined repertoire 
of interactions available for slow control in the central- 
izer, e.g., A = Ug(H) e Zh(G), together with the ones 
capable of supporting fast modulation, e.g., Hamiltonian 
B considered above. Again, the conditions established in 
fll3| provide a necessary and sufficient criterion for uni- 
versality. 

Let us briefly comment on the situation where no 
knowledge is available on the error space £ . In this case, 
decoupling can be achieved only by maximal averaging 
H , so that the effective Hamiltonian is a trivial c-number, 
H e ff — A 1. Since for decoupling groups with this prop- 
erty Z(G) = {1}, control schemes based on multiple or 
twisted decouplers are not useful anymore. In principle, 
one could still attain complete control in two circum- 
stances: either Gb.b. contains a universal set of operations, 
which have to be performed synchronously with the de- 
coupler clock T c ; or a universal set of Hamiltonians can be 
switched fast. Even with this option, the minimum num- 
ber of required b.b. operations, \Q\ = (dim("Hs)) 2 ||, 
may be very large for relevant systems, strength losses in 
effective Hamiltonians (||) becoming possibly quite sig- 
nificant. We turn now to analyze more specifically the 
case of quantum computation. 

Universal Quantum Computation.— Consider a quan- 
tum computer made of K qubits, Hs — (C 2 )® K , and as- 
sume that the relevant coupling to the environment can 
be accounted by a linear interaction of the form 



H. 



SB 



= E^ ; 



B. 



(i) 



x,y,z; i = l,...K , (6) 



for suitable bath operators Ba ■ The above coupling en- 
compasses various models of interest, with error space 
spanned by combinations of singlc-qubit operators. Fol- 
lowing the classification of H], dim(£) = 3K for in- 
dependent decoherence where {E a } — {&a }, while 
dim(£) = 3 in the opposite limit of collective decoherence 
with global generators {E a } = {J2i &i }, intermediate 
situations occurring with cluster decoherence. Selective 
decoupling of the quantum register from H$b requires a 
minimum of 3 b.b. operations i.e., \Q\ = 4. A convenient 
choice is Q — {1, <8>|L 1 a'o }, in which case the g^s corre- 
spond to collective 7r-rotations and Z(Q) is generated by 

bilinear interactions of the form a a Co pjjlq]. 

We sketch now the application of the programming 
schemes described above. Suppose that, in addition 

to the (/-decoupler, we are equipped with a second de- 

(i) 
coupler based on a group Q, where for instance ci 

is interchanged with a z l and a y with ai' (say, i — 
l,j = 2). Obviously, Hg(£) = 0. Then almost any 
set of gates Ua = e lAtA , Ub = e lBtB will be univer- 
sal over H.s, noise-suppression being preserved if Hamil- 
tonians A e Zh{G), B € Zh{G) are applied in paral- 
lel to the decouplers G, G respectively. Since the lat- 
ter groups are connected by a double 7r/2-pulse P = 
exp(— i(-n/A)a y ') exp(— i(ir/4)ai '), the same result can 
be reached if the realizable group of b.b. operations only 
includes Gb.b. = {G, P, P^}, but we have the capability of 
inserting P, P^ pulses synchronized with T c . Finally, in 
the case that full knowledge is available of the decoupling 
sequence, a constructive result can be given. Note that 
the centralizer Zh(G) contains the bilinear Heisenberg 
couplings ^2 a cTaOa enabling one to to implement swap- 
ping between any pair of qubits |lq| . Since controlled- 
not gates can be assembled as a sequence of "square- 
root swaps" and single-qubit operations |16| , universal 
quantum logic can be performed if we have access to fast 
modulation of single-qubit Hamiltonians in addition to 
the required two-qubit interactions in the centralizer. A 
similar conclusion was conjectured in fll7| . 

Whether the proposed approach can be viable in re- 
alistic situations will strongly depend on the details of 
the system and the environmental noise, as well as on 
the sophistication of the available technology. Bang-bang 
control is potentially suitable for NMR quantum compu- 
tation, as long as b.b. pulses are able to be effected with 
a bandwidth small on the scale of the required spectral 
resolutions. At present, an experimental demonstration 
of b.b. control has been reported for all-optical quan- 
tum circuits |l8| |. Despite the challenges involved, the 
appeal of limited space resources may stimulate efforts 
to practically implement decoupling in different quantum 
information processors. 

Conclusion.— We showed how to achieve noise-tolerant 
universal quantum control on the full Hilbert space of 
the system based on purely unitary open-loop manipu- 



lations. From the perspective of quantum information 
processing, this implies the potential of accomplishing 
noise-protected universal quantum computation without 
the cost of extra space resources. The method, which 
is best suited for slow-response non-Markovian quantum 
baths, complements existing approaches based on quan- 
tum error correction, where typically a memoryless error 
scenario is assumed. As a further step toward the goal 
of a truly fault-tolerant computation scheme, the present 
assumption of perfect control resources should be relaxed 
to allow noisy decoupling and programming operations. 
A separate analysis of the issue of robustness shall be 
presented in a future work. 
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